Precise predictions for Higgs production in models with color-octet scalars 
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We describe an effective-theory computation of the next-to-next-to-leading order (NNLO) QCD corrections to 
the gluon-fusion production of a Higgs boson in models with massive color-octet scalars in the (8, l)o representa- 
tion. Numerical results are presented for both the Tevatron and the LHC. The estimated theoretical uncertainty 
is greatly reduced by the inclusion of the NNLO corrections. Color-octet scalars can increase the Standard Model 
rate by more than a factor of two in allowed regions of parameter space. 



1. Introduction 

The hunt for the Higgs boson to uncover its role 
in electroweak symmetry breaking is now being 
undertaken at both the Tevatron and the Large 
Hadron Colhder (LHC). The CDF and DO collab- 
orations at the Tevatron have recently announced 
a 95% exclusion limit on a SM Higgs boson with a 
mass in the range 162 GeV < ruh < 166 GeV [1]. 
This search requires precise theoretical predic- 
tions for the cross section for Higgs production. 
The dominant hadronic production mechanism, 
gluon fusion through a top-quark loop, is known 
exactly through next-to-lcading order in pertur- 
bative QCD j2|3j . In the effective theory with 
TOt — )■ oo, both the NLO corrections (3] and the 
NNLO corrections are known [5 6 .7]. When nor- 
malized to the full mt-dependent leading-order 
result, this effective theory reproduces the exact 
NLO result to better than 1% for rrih < 2mt and 
to 10% or better for Higgs boson masses up to 
1 TeV [5]. The QCD radiative corrections dras- 
tically alter the Higgs production cross section 
prediction; for example, the gluon fusion cross 
section is increased by roughly a factor of three 
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above the LO prediction at the Tevatron after the 
NNLO corrections are included. Only at NNLO is 
an accurate prediction free from debilitating un- 
certainties obtained. After inclusion of the NNLO 
corrections, the remaining uncertainty form un- 
calculated higher-order terms is estimated to be 
roughly ±10%. Updated cross sections for Higgs 
production in gluon fusion for use at the Teva- 
tron and LHC are available in Refs. |9I10| . For a 
recent review of the status of theoretical predic- 
tions for Higgs boson production in the SM, see 
Ref. [11]. 

The properties of the Higgs boson can be mod- 
ified in theories with additional particles, and 
measurement of these properties consequently 
serves as a window into physics beyond the SM. 
As an example of such a Standard Model exten- 
sion, we consider new particles transforming in 
the adjoint representation under the color gauge 
group. Several interesting extensions of the SM 
introduce such states. Scalars in the (8, l)o rep- 
resentation can arise in theories with universal 
extra dimensions |12I13| . The primary decay for 
such states is expected to be into either bb or tt, 
depending on the scalar mass and other model pa- 
rameters. The Tevatron experiments can search 
for these states via pair production of (8, l)o 
scalars, leading to a four 6-jet final state. The 
search reach in the scalar mass was recently es- 
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timated to be 280 GeV 13 . The direct search 
for these scalars is rendered difficult by the large 
QCD background. Reduction of the background 
requires significant cuts that reduce the signal 
and therefore the search reach. It is possible 
that indirect searches for these scalars, such as 
via their infiuence on the Higgs production cross 
section, can probe masses competitive with direct 
searches. 

We summarize here our recent calculation of 
the NNLO corrections to the Higgs boson pro- 
duction cross section in models with a (8, l)o 
scalar [14] . To be as independent as possible from 
the origin of this scalar, we study a simple model 
that couples this scalar to both QCD and the 
Higgs doublet via renormalizable operators. We 
utilize the effective-theory approach valid when 
both the SM top quark and the new scalar are 
heavier than the Higgs boson. As a byprod- 
uct of our computation, we derive all renormal- 
ization constants needed to study higher-order 
corrections in models with these scalars. We 
use the renormalization-group equations govern- 
ing the evolution of the scalar-sector couplings to 
determine the likely ranges of the various param- 
eters which appear. 

We also study the phenomenological implica- 
tions of the color-octet scalar for the Higgs pro- 
duction cross section at both the Tevatron and 
the LHC. Only at NNLO is the scale dependence 
sufficiently reduced to allow precise predictions 
for the scalar-induced effects to be obtained. The 
deviations from the Standard Model prediction 
for Higgs production are large at both collid- 
ers. Deviations of a factor of two, much larger 
than the estimated uncertainties from higher- 
order corrections and parton distribution func- 
tion errors, are obtained for scalar masses near 
the estimated direct search reach of 300 
GeV. 

2. Details of the Model 

We begin with the following Lagrangian, which 
describes the Standard Model coupled to a color- 
octet scalar in the (8, l)o representation: 



-gl G4S Tr [S^] ^ - XiH^H Tr [S^] 
-A,(^i7t^-^)'. (1) 

S denotes the matrix-valued scalar field S = 
gArpA^ indicates the Higgs doublet before 
electroweak symmetry breaking, v is the Higgs 
vacuum-expectation value, and I?^ is the covari- 
ant derivative for adjoint fields. The Higgs quar- 
tic coupling has been explicitly included to de- 
fine its normalization. After electroweak sym- 
metry breaking, the Higgs doublet is expanded 
as H = (^0,{v + h)/V^ in the unitary gauge. 
The mass of the color-octet scalar becomes m| = 
rrig + Aiu^/2. The Feynman rules which describe 
the scalar couplings to the Higgs boson h and 
to gluons are easily obtained from Eq. ([T|). The 
free parameters which govern the scalar proper- 
ties are ms, Ai, and G45. We note that higher- 
order operators which break the — > — S" sym- 
metry present above and which allow the scalar 
to decay are obtained in explicit models which 
contain this state |12ll3j . We neglect them here 
since we anticipate that they have little effect on 
the gg h production cross section. We com- 
ment on the appearance of the Tr [S^] term m 
the scalar potential. A quartic-scalar coupling is 
generated by QCD interactions even if it is set 
to zero at tree-level. At NNLO the quartic cou- 
pling must be included to obtain a renormaliz- 
able result. We include this operator in the tree- 
level Lagrangian with a coefficient scaled by g^, 
the QCD coupling constant squared, to permit an 
easier power-counting of loops. 

3. Calculational details and analytic re- 
sults 

We sketch here our derivation of the effective 
Lagrangian describing the interaction of the Higgs 
boson with gluons through NNLO in the QCD 
coupling constant. When both the adjoint scalar 
and the top quark are heavier than the Higgs bo- 
son, they can be integrated out to derive the fol- 
lowing effective Lagrangian: 



Z:^"" = CsM +Tr[D^SD''S]-m'^Tr[S'^] 



(2) 
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where Ci is a Wilson coefficient and the operator 
d is 

1 

4 

The effective Lagrangian for the gluons and light 
quarks, Cqq[/ , has the same form as Cq^j^ ex- 
cept that the fields and parameters it contains 
are rescaled by decoupling constants that account 
for the effects of the heavy states. In order to 
distinguish the fields and parameters occurring 
in the full Lagrangian from those in the effec- 
tive Lagrangian, we denote the latter ones with 
a prime. In this manuscript we integrate out 
the top quark and the color-octet scalar in a sin- 
gle step. Our calculation is therefore a two-scale 
problem. This is clearly demonstrated by the ex- 
ample diagrams contributing to the NNLO Wil- 
son coefficient shown in Fig. [T] 



s >- 




Figure 1. Example two-scale diagrams contribut- 
ing to the Wilson coefficient at NNLO. 



The Wilson coefficient Ci can be constructed 
by computing the amplitude for gg ^ h in the 
limit that the initial gluon momenta vanish, as 
reviewed in Ref. fTSl. We start from the relation 



V {N-^-l){d~2){p^-p2Y 

X A^°;rMPi,P2)|p,=p.=o (4) 

between the bare Wilson coefficient and the 
bare amplitude for the process gg ^ H in the 



full theory. Here, pi and p2 are the momenta 
of the two gluons, N is the number of colours, 
and d = 4 — 2e is the dimension of space-time. 
The factor is the bare decoupling coefficient by 
which the bare gluon fields in the full and effective 
theories are related: 



G' 



^0 ^0,a 



(5) 



It can be expressed as 

C3" = i + n^(p-o), (6) 

where Ilg. is the transverse part of the bare gluon 
self-energy in the full theory. In order to obtain 
the Wilson coefficient through NNLO in the QCD 
coupling constant, M.^^'^^^^ is needed up to three 
loops, while Ii%{p = 0) is needed through two 
loops. We can set all scaleless integrals to zero in 
dimensional regularization, and as a result, only 
diagrams containing at least one massive scalar 
propagator contribute to ^^q{p — 0) and to the 
right hand side of Eq. ^ . We perform our calcu- 
lations in a covariant gauge with gauge parameter 
^, leading to the following gluon propagator: 



(7) 
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In all diagrams, terms up to first order in ^ are 
retained. All ^-dependent terms cancel in TIq{p = 
0) up to two-loop order, and in A^^^'^i^^ up to 
three-loop order, providing a strong calculational 
check on the result for Cj*. 

We now outline the steps involved in deriving 
the Wilson coefficient describing the color-octet 
interactions with the Higgs boson. Initially, all 
quantities appearing in Eq. Q are expressed in 
terms of the bare masses m^, m^, and the bare 
coupling constants g'^, G^g. We renormalize the 
scalae and top-quark masses in the MS scheme 
according to 

"^s = V^ms "^s, rn'l = m^. (8) 

The explicit expressions for these renormaliza- 
tion constants, and all others, can be found in 
Ref. [14j . The quartic coupling first appears in 
the gg ^ h amplitude at the two-loop level, and 
its renormalization is therefore only required to 
one loop: 

G% = Zis G4S ■ (9) 
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We convert the bare strong coupling constant of 
the full theory into the bare coupling of the effec- 
tive theory using decoupling constants obtained 
from the ghost self energy and the ghost-gluon 
vertex up to two-loops, as described in detail 
in [16115] . We then renormalize it in the effec- 
tive theory using the MS scheme: 

gf ^^,^Z'^g',. (10) 

The Wilson coefficient itself requires a renormal- 
ization factor [17118119] : 

Ci = ^C',. (11) 
'ill 

Our final result for the renormalized Wilson co- 
efficient is rather lengthy, and can be found in 
Ref. [M]. 

4. Numerical Results 

We now present numerical results for the Higgs 
boson production cross section in gluon fusion, to 
study the deviations induced by the color-octet 
scalar. We include the effects of the top quark, 
the scalar, and also the bottom quark on the 
Higgs cross section. We comment first on pre- 
cisely what terms we include in the cross sec- 
tion. We denote by gt+s the terms obtained 
by squaring together the top and scalar ampli- 
tudes, and keeping both the interference term and 
the pieces from each separate particle squared. 
We let aTBjO'sB denote the interferences between 
the bottom-quark amplitude with the top and 
the scalar pieces, respectively. asB indicates the 
contribution from the bottom-quark amplitude 
squared. For the cross section at the 7i-th or- 
der in perturbation theory, we use the following 
expression: 

cr" = aT'!ls{mt,ms) K^FT + <^SBi^S,mb) 
+ a^%{mb) + cJ^%{mt,mb). (12) 

K^prp denotes the ratio of the n-th order cross 
section over the LO result, with both computed in 
the effective theory defined in Eq. ([2]). The cross 
section multiplying the if-factor is the LO cross 
section with the exact dependence on the scalar 
and top-quark masses. The remaining terms 



account for the scalar-bottom interference, the 
top-bottom interference, and the bottom-squared 
contribution at LO with their exact mass depen- 
dences. In the SM, the scaling of the exact LO 
cross section by the EFT X-factor is known to 
furnish an approximation accurate to the few- 
percent level or better for Higgs masses below 
roughly 400 GeV For the scalar, this ap- 
proximation has been studied against the exact 
NLO calculation [20], and is again accurate to 
the 1 — 2% percent level for rrih < ms, which 
is the region we focus on here. If desired, the 
exact NLO corrections to the top-bottom inter- 
ference and bottom-squared terms 3 21j can be 
included, as can those to the bottom-scalar in- 
terference [20]. These affect the cross section at 
the 1 — 2% level, and for simplicity are neglected. 
Various electroweak corrections known for the SM 
contribution [22|23|9|24] are not known for the 
scalar, and for consistency are neglected. 

We use the pole mass nit — 173.1 GeV for the 
top quark [25] , while for the bottom quark we 
use the MS mass [26| • The choice of pole or AIS 
mass for the 5-quark has an insignificant effect 
on the fractional deviation between the scalar- 
induced cross section and the SM result. We 
use the MSTW parton distribution functions [27] 
extracted to the appropriate order in perturba- 
tion theory. For the scalar sector, we must set 
the parameters Ai, G4S, and tos- We perform 
a renormalization-group analysis to constrain the 
possible values of the scalar couplings. We show 
in Fig. [2] the maximum allowed value of Xi(v) 
obtained by demanding Ai(/i) remain perturba- 
tive until n ~ 10 TcV, as a function of the Higgs 
mass, for the choices G^siv) — 1 and G4s{v) = 0. 
For a complete description of the renormalization- 
group analysis, we refer the reader to Ref. [14] . 

5. Results for the Tevatron and the LHC 

We now present numerical results for both the 
Tevatron with = 1.96 TeV and the LHC, for 
which we assume ^/s = 7 TeV. We begin by show- 
ing in Fig. [3]the LHC cross section for ms = 300 
GeV at LO, NLO and NNLO in QCD perturba- 
tion theory, to see the effect of including higher- 
order QCD corrections. The renormalization and 
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Figure 2. Maximum value of Ai allowed by 
perturbativity for the choices G4s(?;) = 1 and 
G4s{v) = 0, as a function of Higgs mass. 



factorization scales are equated to hf = IJ'R = IJ', 
and are varied in the range ruh/i < fj. < mh, con- 
sistent with previous studies of the Higgs produc- 
tion cross section in the SM. The scale-variation 
errors are large at both LO and NLO, and the 
corresponding error bands do not overlap. Only 
at NNLO can a reliable prediction for the cross 
section be made. It is also clear from this fig- 
ure that large variations from the SM prediction 
are possible for scalar masses near the expected 
Tevatron limit of mg w 300 GeV. The cross sec- 
tion differs from the SM result by more than a 
factor of two for this parameter value. 

To study further the effect of the color-octet 
scalar on the Higgs cross-section prediction at 
the Tevatron, we show below in Fig. 0] the cross 
sections for the example Higgs mass ruh = 165 
GeV as functions of the scalar mass. Deviations 
from the SM are visible over scale errors for scalar 
masses of up to 1 TeV. The scalar contributions 
to the Higgs production cross section are large, 
and direct searches are hindered by the need to 
pair produce the scalars and by the large QCD 
background. This suggests that the indirect con- 
straints on the scalar parameter space coming 
from the Tevatron Higgs exclusion limit could be 
as strong as the direct search reach. 



Figure 3. Higgs production cross section in gluon- 
fusion at the LHC for ms — 300 GeV as a func- 
tion of the Higgs boson mass. The bands indi- 
cate the scale variation mh/A < < mh- From 
bottom to top, the bands indicate the variations 
of the LO, NLO, and NNLO cross sections. Ah 
other parameters are as described in the text. For 
orientation, the SM result at NNLO for the cen- 
tral value /i — mfi/2 is shown. We note that the 
residual scale uncertainty of the SM cross section 
is approximately ±10%. 



6. Conclusions 

We have discussed a study of the effects of 
(8,l)o scalars on the Higgs-gluon effective La- 
grangian. We have presented phenomenological 
predictions for both the Tevatron and LHC. The 
scalar-induced effects are large at both collid- 
ers, and are visible over both scale and PDF er- 
rors for scalar masses of up to 1 TeV. This sug- 
gests that the scalar parameters can be strin- 
gently constrained using the exclusion limit on 
the Higgs boson production cross section ob- 
tained by the Tevatron, and that the limits would 
be competitive with those obtained from direct 
searches for scalar pair production. Models that 
include color-octet scalars are an interesting and 
phenomenologically-rich example of physics be- 
yond the SM. The NNLO calculation presented 
here quantifies precisely the effect of a (8, l)o 
scalar on the Higgs production cross section in 
gluon fusion, and our phenomenological study 
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Figure 4. Higgs production cross section in gluon- 
fusion at the Tevatron for nih — 165 GeV as a 
function of the scalar mass. The bands indicate 
the scale variation mh/4: < fj. < nih- All other 
parameters are as described in the text. Also 
shown is the SM cross section with its correspond- 
ing scale uncertainty. 



shows that the deviations induced by the scalar 
are large at the Tevatron and the LHC. The ex- 
clusion limit on Higgs production set by the Teva- 
tron collaborations allows for the scalar parame- 
ters to be stringently constrained. 
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